Search for Small-Mass Black Hole Dark Matter with Space-Based Gravitational Wave 

Detectors 



O 

o 



> 

(N ■ 

vn ■ 

o : 
o ■ 

Or 

6 ■ 

a: 

o3 ■ 



Naoki Seto and Asantha Cooray 
Theoretical Astrophysics, Division of Physics, Mathematics and Astronomy, 
California Institute of Technology, MC 130-33, Pasadena, CA 91125; (seto,asante)@caltech.edu 

The high sensitivity of upcoming space-based gravitational wave detectors suggests the possibility 
that if halo dark matter were composed of primordial black holes (PBHs) with mass between 10 16 g 
and 10 20 g, the gravitational interaction with detector test masses will lead to a detectable pulse-like 
signal during the fly-by. For an improved version of the Laser Interferometer Space Antenna with 
a reduced acceleration noise at the low-end of its frequency spectrum, we find an event rate, with 
signal-to-noise ratios greater than 5, of ~ a few per decade involving black holes of mass ~ 10 17 g. 
The detection rate improves significantly for second generation space based interferometers that are 
currently envisioned, though these events must be distinguished from those involving perturbations 
due to near-Earth asteroids. While the presence of primordial black holes below a mass of ~ 10 16 g 
is now constrained based on the radiation released during their evaporation, the gravitational wave 
detectors will extend the study of PBHs to a several orders of magnitude higher masses. 



Based on observed rotational velocity measurements of 
the Milky Way disk, the presence of a substantial dark 
matter component in the Milky Way halo is now well 
established M. With no unique guidance as to the na- 
ture of dark matter, a large number of candidates based 
on both astronomical arguments, such as baryonic dark 
matter involving substellar mass remnants , and parti- 
cle physics expectations are now routinely considered 
to explain the missing mass. An interesting possibility 
is that the halo dark matter is composed of primordial 
black holes (PBHs) 4]. While PBHs are expected over 
a rather wide range of masses 0, the population below 
~ 5x 10 14 g is expected to have evaporated by the Hawk- 
ing radiation over the age of the universe , while holes 
with a mass slightly above this limit are still emitting 
high energy particles today 0. With observations re- 
lated to the galactic Gamma-ray background, the evapo- 
rating black hole contribution to the dark matter density 
in the Universe, in terms of the critical density, is at the 
level roue hly below 1(T 8 0. 

From the high mass end, constraints on the possibil- 
ity that dark matter is composed of PBHs come from 
dynamical arguments such as the potential disruption of 
galactic clusters and similar bound structures |8j . These 
constraints generally limit the primordial black hole mass 
to be roughly below few solar mass. Similar constraints 
come from galactic microlensing observations, which lim- 
its primordial black hole masses to be below roughly few 
tenths solar masses @. In combination, at least over 
a 10 16 decade in primordial black hole mass, between 
evaporating holes with mass below 10 17 g and ~ 0.1 M Q , 
remains yet to be studied and techniques to understand 
their presence are limited (see also 10]). Though the 
galactic gamma-ray background and its weak anisotropy 
suggests a evaporating PBH density of roughly 10 10 pc~ 3 , 
PBHs can explain the total halo dark matter if the PBH 
mass spectrum enhances the abundance above 10 17 g [Tlj . 
A promising possibility of probing small-mass PBHs in- 
volves femtolensing of gamma-ray bursts (GRBs) both 



in terms of an interference pattern in the frequency spec- 
trum |12| or relative flux differences in Astronomical Unit 
( AU)-scale separated light curves [l3( . The existing data 
weakly constrain the mass density fraction, relative to the 
critical density, of compact objects with mass between 
10 17 g and 10 20 g to be below 0.2, if the average redshift 
of GRBs is unity [14j . This technique only constrain the 
extragalactic density, though, it would be interesting to 
see if the prediction [Tl| related to the galactic dark mat- 
ter can directly be tested. 

The PBHs, especially at the small-mass end, are not 
expected to interact or be captured by other massive bod- 
ies such as the Sun; they can be considered as another 
weakly interacting massive particle (WIMPs). In con- 
trast to active PBHs with masses < 10 15 g, the physical 
nature of PBHs with > 10 17 g is characterized only by 
their masses and extremely small size. In reality, they 
can be detected only through the gravitational interac- 
tion. But, as their masses are much larger than standard 
WIMPs predicted by particle physics, the expected flux 
of PBHs would be very small and we need detectors with 
a large effective area for their search. 

In this Letter, we suggest that PBHs around a mass 
of 10 17 g could in fact be detected directly with planned 
space-based gravitational wave detectors that have large 
cross sections to moving massive bodies. These missions 
include the Laser Interferometer Space Antenna (LISA) 
flEf and several vision missions of future that are now 
routinely discussed in the long term planning by var- 
ious space agencies: Gravitational Echos Across Time 
(GREAT; pjf ), the Big Bang Observer (BBO) of NASA's 
Structure and Evolution of the Universe Roadmap, and 
a similar mission concept named Decihertz Interferome- 
ter Gravitational Wave Observatory (DECIGO; [T3) (see 
also El)- The direct detection simply involves the 
gravitational interaction between the fly-by PBHs and 
detector test masses such that with high sensitivity even 
a small gravitational perturbation can eventually be seen 
above the detector noise. As an example, the LISA mis- 
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sion is expected to probe the frequency range between 
1CP 5 Hz and 10 _1 Hz and detect perturbations above 
h ~ SL/L ~ 10 ~ 21 at ~ lmHz with detectors separated 
by arm lengths L ~ 5 x 10 6 km. 

The interaction of a PBH and a gravitational- wave ob- 
servatory test mass can be considered either as a direct 
interaction, when the fly-by separation is smaller than 
the arm length, or as a tidal interaction, in the case where 
arm length is smaller than the fly-by distance. We first 
consider the former situation and take the scenario with 
a PBH mass M, passing by one of the test masses of the 
interferometer with a relative velocity V with the clos- 
est approach distance R. The time dependent accelera- 
tion a(t) of the spacecraft test mass towards the direction 
of the closest approach has a single pulse-like structure 
given by 



a(t) = 



GM.R 



[R 2 + (Vt) 2 } 3 / 2 ' 



(1) 



where we have set the origin of the time coordinate, t = 0, 
to be at the closest approach. Note that the fly-by PBH 
can also perturb other test masses of the interferome- 
ter configuration since at least 3 or more spacecrafts are 
generally utilized in future gravitational wave observato- 
ries. When making an order-of-magnitude estimate we 
neglect additional information related to such perturba- 
tions and ignore aspects that involve the relative con- 
figuration of interferometer-body system related to the 
perturbing PBH. 

The Fourier component a(/) of eq. is given in terms 
of the 1st order modified Bessel function K\ as 



/DO (~t yr 
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For the direct detection, the signal-to-noise ratio (SNR) 
of the pulse a(t) is formally written as SNR 2 = 
2 Jo °(/) 2 / a n(f) 2 df with the noise spectrum a n (f) of 
the detector. In the case of space-based interferometers, 
the noise has two main contributions involving the proof- 
mass noise and the optical path noise with the former 
dominating the low frequency regime. As we will soon 
discuss, the typical signal a(f) due to the fly-by PBH of 
mass around M, > 10 17 g has support in the SNR inte- 
gral when / < 10 -4 Hz. This is in the the low frequency 
regime of all the proposed future interferometers, includ- 
ing the LISA mission. The proof-mass noise of LISA is 
estimated to be constant down to ~ 10~ 4 Hz; It might 
become a factor of two larger at 3 x 10~ 5 Hz and more 
at lower frequencies |15| . however, the low frequency be- 
havior of the noise curve, around 10~ 5 Hz, is yet to be 
determined precisely. In order to estimate the SNR, we 
simply extrapolate the constant proof-mass noise down 
to / = 0, and set a n {f) — 2a = const, where a is the 
proof-mass noise for a single inertial sensor in various 
interferometers. While frequency dependence of noise 
below 10 -4 Hz is ignored, we do pay attention to the 
frequency dependence of the fly-by pulse event so as to 



understand requirements on the low frequency behavior 
of the proof-mass noise. 

After some straight forward algebra, we obtain the 
SNR as 



SNR 2 = 



3tt (GM.f 
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Note that this relation corresponds to the second for- 
mula of eq. (4.35) in Ref. 0], but our result differs by 
1/13 ~ 8% due to the exact use of a limit associated 
with the Bessel function, instead of an approximation 
used there, and after taking in to account a difference in 
the definition related to the proof-mass noise. 

Since the above expression is for the case where 
the closest approach distance R to a single spacecraft 
is smaller than the arm-length L of the interferometer, 
namely R < L, we now consider the opposite case with 
R > L. Now, the tidal deformation of the interferome- 
ter is a measurable effect, but this involves an additional 
suppression factor ~ n(L/R) in eq. l(T|. where we have 
ignored the numerical prefactor n as we do not consider 
detailed geometry of these events, though we note that 
this factor could be greater than unity (up to 2) in cer- 
tain favorable configurations of the interferometer with 
respect the trajectory of the PBH fly-by. In this tidal 
limit, we find the SNR to be 



SNR 2 - 3 * (GM -^ 
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This is the relevant expression for LISA and 
BBO/DECIGO in the interesting mass range 
M. > 10 17 g, though, at the low mass end (10 16 
g), the former expression applies for LISA. 

Using eq. J3J we can now estimate the maximum length 
of the closest approach i? max for a given SNR threshold 
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where ciq — 3 x 10~ 15 m/s 2 \/Hz and L = 5 x 10 11 cm 
are the reference fiducial parameters for LISA in its cur- 
rent design 0] (at least the proof-mass noise around 
/ ~ 10~ 4 Hz). We have taken the typical velocity disper- 
sion V = 220-\/5/2 = 350 km/s of the halo dark matter 
particles relative to the solar system using the Galactic 
rotation velocity of 220 km/s and the Galactic radius to 
the Solar system of r g ~ 8 kpc following Ref. @. The 
gravitational perturbation involves a pulse like signal in 
the data streams with a characteristic frequency given by 
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As mentioned already, for most planned space interferom- 
eter such as LISA, such a low frequency is in the regime 
where the detection is limited by the proof-mass noise. 

When we assume that the halo dark matter with den- 
sity pdm ~ O.OIIM0/PC 3 [2l| around the Sun is made 
with PBHs of mass M«, their flux, F — pdmVM" 1 , is 

F = 9.0 x 10 27 (o o^Qpc-s) (350km/s) 

and the fly- by event rate, rj — ttFR^^, above a certain 
SNRis 

• nm/ M ' Y 1/5 ( SNR\- 4/5 ( V \ 3/5 

x / PDM \ (±Y A/5 (L_\ i/b (N\ _! 

X ^O.OllM pc-V \a J \Lo) \ 1 J Y * ' 

(8) 

where N represents the effective number of interferom- 
eters with sufficient relative distances. While LISA in- 
volves a single set of interferometers, some planned mis- 
sions, e.g. BBO, plans to use multiple interferometer 
arrays with large separation to improve the localization 
of binaries sources such that N > 2. If the signal is to be 
extracted using a correlation analysis, especially for low 
signal-to-noise events, of two arrays, there is no further 
increase in the event rate and one should continue to use 
N = 1. 

For the reference parameters of LISA, assuming a use- 
ful SNR threshold of 5, the detection rate is ~ 0.1 per 
ten years as shown in eq. iJSJ. In the tidal limit, the 

— 1/5 

event rate weakly depends on mass M, by rj oc M % 
The combination (aL^ 1 ) that appears in all of the above 
expressions determines the low frequency sensitivity of 
the interferometers to the gravitational wave amplitude 
h. Therefore, once the threshold distances i? max are con- 
firmed to be larger than the arm- length L, we can eas- 
ily compare the event rate for various interferometers by 
studying their noise curves at the low frequency regime. 
In the case of LISA, when the mass M. falls below 10 17 
g, the maximum fly-by distance approaches that of the 
arm length (~ 5 x I0 11 cm). In this case, we repeat the 
event rate using the SNR involving the direct perturba- 
tion (eq|3J); in this limit, for reference, the event rate 

1 /3 

depends on the mass M. by 1) oc M m ' . 

In figure I, we show the potential detectability of 
fly-by events involving PBHs with mass M, using var- 
ious planned space-based interferometers. Instead of the 
event rate, our results are presented in terms of the local 
halo density of PBHs required for a detection, with SNR 
greater than 5, of one event per decade. For reference, to 
put these detections in the context of the galactic dark 
matter density, we also show current observational con- 
straints for the PBH density from the 7-ray background 
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FIG. 1: The detection thresholds for PBH fly-by with various 
upcoming space-based gravitational wave observatories (in 
solid lines from top to bottom: LISA, GREAT intermediate- 
frequency-mission, BBO/DECIGO, and GREAT very-low- 
frequency-mission). The thick solid line shows the expected 
density of halo dark matter in the form of black holes; at 
M, — 10 15 g, the density is constrained by 7-ray background 
observations such that the particle density is below ~ I0 10 
pc -3 , while no similar constraints exist above M. •> I0 17 
g and we take the whole halo to be formed of black holes 
with the mass given in the horizontal axis (with a density 
O.Off Mq pc~ 3 ). The detection limits assume an event rate 
of a single detection per decade; if no events are detected, 
these curves would roughly correspond to the constraint one 
can put on the black hole contribution to the halo dark mat- 
ter density. In each of the instruments considered, we la- 
bel the following three parameters: (arm-length/proof-mass 
noise/total number of detector arrays) relative to those refer- 
ence values of LISA (arm length: 5 x 10 6 km, proof-mass noise: 
3 x 10 _15 m/s 2 /\/Hz). The symbols on each of the curves rep- 
resent the average frequency of the gravitational perturbation 
produced on the interferometer by a black hole of correspond- 
ing mass with triangle, diamond, circle and square represent- 
ing frequencies of 10 -7 ,10 -6 ,10 -5 and 10 -4 Hz, respectively. 
In the case of LISA, the threshold distance R equals to the 
arm length L at M. = 0.87 x 10 17 g where the slope of the 
curve changes. 



(M, < I0 17 g) an d the expected density of dark mat- 
ter based on a Galactic model (M. > 10 17 g) where the 
density is now determined to be (0.011 ±0.005) M Q pc~ 3 
[2l|. In the case of LISA, the transition (R = L) between 
eqs. (3) and (4) occurs at M. = 0.87 x 10 17 g, but for all 
other interferometers, we find R > L in the mass range 
considered in figure 1. 

In figure 1, next to the name of each of these in- 
terferometers, we also denote three basic numbers that 
characterize an interferometer (arm-length/proof-mass 
noise/number of detector arrays) normalized by the refer- 
ence values for LISA. The GREAT-intermediate mission 
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|l|j has a marginal possibility for a PBH detection at 
M, ~ 10 17 g while The GREAT- very low frequency mis- 
sion with parameters (50/4 x 10~ 4 /1) would detect them 
even if less than a few percent of the halo dark matter 
is made with PBHs with mass between 10 17 g and 10 20 
g. Incidently, the low-frequency end of the LISA noise 
has not been finalized and there are strong scientific rea- 
sons to improve it beyond what was originally suggested 
|20j . If we take the most improved low-frequency noise 
from Ref. |2(j with scaled parameters (5/0.1/1), we find 
a detection level similar to that shown for the GREAT- 
intermediate mission suggesting that even with an im- 
proved version of LISA, a detection with an event rate of 
few per decade, at mass of M, ~ 10 17 g, may be within 
the reach of first-generation interferometers. We note 
here that if the event rate is too high with > 100 yr -1 , 
as would be the case for GREAT-low frequency mission 
|16|. the signals will be overlapped in frequency space 
due to the low frequency nature. This will complicate 
the extraction and would involve a detailed analysis sim- 
ilar to the ones proposed to extract individual gravita- 
tional waves in their confusion limit. In this figure, we 
also note the characteristic frequency of the fly- by event. 
This is done to emphasize the very low frequency aspect 
of the detection, especially with respect to the noise re- 
quirements associated with proof-mass, or acceleration 
noise. 

In the search for fly- by PBHs, it is important to distin- 
guish their signals from other signals in the interferom- 
eter. The primary targets of a space-based interferom- 
eter are low frequency gravitational waves. In general, 
their measurement is limited to some extent by the as- 
trophysical confusion noise, depending on the frequency. 
For example, the LISA detector noise is expected to be 
dominated by the Galactic white dwarf binary confusion 
noise in the frequency band between 0.1 mHz < / < 3 
mHz |l9j . At a lower frequency than this range, gravi- 
tational waves from merging massive black hole binaries 
might form a confusion noise |19|| . When we measure 
the local acceleration related to PBH fly-by events, how- 
ever, we can largely reduce the relative contribution of 
the low frequency gravitational waves (with wave-length 
^gw > L) in the data streams by taking a certain com- 
bination of the data stream now knowns as the Sagnac 
|22j . This data combination nearly cancels gravitational 
waves using the fact that they propagate with speed of 
light, though it would still be affected by fly-by pertur- 
bations. Therefore, the separation of PBH signals with 
other gravitational wave signals is not as severe as the 
confusion noise problem for gravitational wave observa- 
tions alone. 

There is one source of confusion, however. This in- 
volves similar pulses produced during the passage of near- 
Earth asteroids close to detectors (see, section 4.2.5 of 
Given the expected flux of minor bodies in the so- 
lar system, as determined by various observational data, 



the asteroid perturbations are dominated by those at the 
high mass end between 10 13 g to 10 15 g, or sizes around 
a km or slightly less, with an event rate of the order 
0.05 yr~* (the rate for smaller mass events are substan- 
tially less and can be ignored). Unfortunately, though 
the relative velocity of asteroid events are smaller (~ 30 
km sec -1 ) than those involving PBHs, these events could 
have similar frequencies (as PBHs) due to differences in 
the maximum distance to which they can be detected. 
On the positive side, orbital parameters of roughly 10% 
of such near-Earth asteroids are already known while this 
fraction is soon expected to grow substantially with dedi- 
cated near-Earth asteroid search programs [23j such that 
during the operation of these gravitational wave obser- 
vatories, expected perturbations from asteroids can be a 
priori determined. Moreover, we expect optical followup 
searches of all candidate events to further refine the sam- 
ple and to confirm the compactness of perturbers. Detec- 
tion of a fly-by pulse with no optical counterpart would 
be a minimum criteria for the selection of a PBH event. 
For this perpose, we need a systematic optical survey 
around the interferometric area above some magnitude 
threshold. With increasing sensitivity, such as the case 
with GREAT-low frequency mission, event rate of aster- 
oids increases such that they might become the primary 
confusion noise; Removal of this confusion noise requires 
independent estimates of the asteroid mass at the few 
percent level in addition to the orbit. 

While with gravitational wave detectors the presence 
of PBHs can be established, it is not easy to determine 
masses of individual events from a single pulse signal that 
is characterized mainly by two numbers: the amplitude, 
M/R 3 , and the time scale, R/V, made from three variales 
M, R and V. This is due to the fact that there is an 
unknown associated with the fly-by distance or the PBH 
trajectory. If multiple detectors are perturbed, one can 
establish the trajectory and then used that information 
to determine mass. 

While there is some possibility for PBH studies with 
LISA, the possibilities are greater with second generation 
space-based missions for which concept studies have al- 
ready begun. For example, BBO is designed to detect in- 
flationary generated gravitational wave background with 
amplitude Slew 1^ 10~ 16 around 1 Hz where the white 
dwarf binary confusion noise is absent. The relevant fre- 
quency for PBH search (/ < 10 _5 Hz) is largely different 
from 1 Hz. However, if the mission design is such that the 
proof-mass noise can be controlled well below ~ 10 _5 Hz, 
such a mission would become a powerful tool for PBH 
search. 
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